Minimal Synthesis of Two-Variable 
Reactance Matrices* 
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(Manuscript received July 3, 1968) 

A simple algebraic method stemming from ideas in minimal state-variable 
realization theory is developed for the synthesis of two-variable reactance 
matrices. The method rests mainly on the factorization of a one variable 
polynomial matrix which is para-Hermitian and positive semidefinite on 
the imaginary axis, and always yields a realization minimal in both vari- 
ables. 

I. INTRODUCTION 

Two-variable reactance functions and matrices, originally intro- 
duced to represent the characteristics of lumped passive networks 
with variable elements, 1 - 2 have become more important because of 
their application to the synthesis of lumped-distributed networks. 
Ansell first showed the two-variable reactance property of networks 
composed of lossless transmission lines and lumped reactances. 3 The 
two-variable theory has also been applied to the synthesis of net- 
works consisting of lumped resistors capacitors and uniformly dis- 
tributed RC lines, 4 - 5 which are of importance in microelectronic struc- 
tures. 6, 7 Besides the various applications, the two-variable reactance 
theory is of theoretical interest in itself since it can be shown that 
passive RLC synthesis is a special case of two-variable reactance 

synthesis. 2 

Koga 8 demonstrated that every nXn two-variable reactance matrix 
W {p, s) can be realized as the impedance seen at the first n ports of 
a lossless (n+qr) -port network in the p-plane terminated at its last 
qr ports with unit inductors in the s-plane; q is the rank of W{p, s), 
and r is the highest degree of s in the least common denominator of 



* This work is based on Chapter III of the author's dissertation, "Synthesis of 
Lumped-Distributed RC Networks" submitted in partial fulfilment of the re- 
quirements for the Ph.D. degree at Stanford University, May 1967. 
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the elements of w. The method is quite complicated and rests heavily 
on the theory of algebraic functions and the structure of para unitary 
matrices. Also it does not guarantee the use of a minimum number of 
elements. Youla 9 solved the problem of synthesizing a lossless two- 
variable scattering matrix by adapting an earlier method for syn- 
thesizing one-variable scattering matrices. 10 The method could be 
adapted to the direct synthesis of an impedance matrix but appears 
to be unduly complicated because of the need to find the transforma- 
tion required to transform a generally unrealizable coupling network 
into a realizable one. 

A simple algebraic method stemming from ideas in minimal state- 
variable realization theory and having similar beginnings as that of 
Youla 9 is developed here for the synthesis of two-variable reactance 
matrices. The method rests mainly on the factorization of a one- 
variable polynominal matrix which is para Hermitian* and positive 
semidefinite on the imaginary axis. Such a factorization is well known 
in n-port network theory and once it is accomplished, the coupling 
network is obtained by simple matrix operations. Furthermore the 
method always yields a network minimal in both types of elements. 

We first introduce some basic definitions and necessary theorems, 
and later we add more as the need arises. The synthesis procedure is 
developed in Section III. Since the various proofs involved are rather 
indirect and tend to cloud the simplicity of the actual procedure, the 
synthesis procedure is outlined in Section IV. The reader interested 
only in the procedure and not in the theory behind it may go directly 
to Section IV where step-by-step instructions are given for the syn- 
thesis of any two variable reactance matrix. In Section V an example 
is worked out. The notation used in this paper is almost the same as 
found in earlier work to assure easy reading for those familiar with it. 9 
Capital letters indicate matrices; bold face letters indicate matrix 
transposition. A superscript dagger indicates the substitution of — s 
or — p for s and p respectively, in the case of two-variable functions. 

II. BASIC DEFINITIONS AND THEOREMS 

The basic notion in the two variable theory is that of a two variable 

positive real matrix, which is a straightforward extension of the same 

notion in the one variable theory (See p. 96 of Ref. 11 and p. 32 of 

Ref. 8). 

* A matrix A(p) is said to be para-Hermitian if A(p) = A*(p) where the bold 
face letter denotes matrix transposition and the superscript dagger denotes re- 
placement of p by —p. 
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Definition 1: An n X n matrix W(p, s) is said to be a two variable 
positive real matrix if 

(i) W is real for real p and s. 
(ii) W is analytic in the domain Re p > and Re s > 0. 

(Hi) W + W* is positive semidefinite in the domain Re p > and 
Re s > 0. 

By statements such as: "W is analytic" in the definition and in 
what follows, we mean, "each element of W is analytic." A two vari- 
able function is said to be analytic at a point if it has a total differ- 
ential at the point. The bold face letter indicates matrix transposition, 
and the superscript star indicates the complex conjugation of each 
element. 

If W{p, s) satisfies conditions (ii) and (iii) of the definition and not 
necessarily condition (i), it will be called a two variable positive 
matrix. 

Definition 2: An n X n matrix W(p, s) is said to be a two variable 
reactance matrix if 

(i) W is a two variable positive real matrix. 

(ii) W + W* = 0. 

The superscript dagger indicates the operation of substituting — p 
and -s for p and s in the original matrix. This definition of a two 
variable reactance matrix is similar to the corresponding one in the 
one variable theory. (See p. 102 of Ref. 11 and p. 32 of Ref. 8.) Anal- 
ogously, as in the one variable case (p. 117 of Ref. 11), it is generally 
hard to check if condition (iii) of Definition 1, which involves the 
whole domain Re p > and Re s > 0, is satisfied for a given two 
variable matrix; we would like to find an equivalent set of conditions 
that are easier to check. In the case of two variable reactance matrices, 
the following theorem proved by Ozaki and Kasami 2 in the scalar 
case, and extended to nonsymmetric matrices by Koga, (p. 33 of Ref. 
8) serves this purpose. 

Theorem 1: The necessary and sufficient conditions for ann X n matrix 
W(p, s) to be a two variable reactance matrix are: 

(i) W is rational in p and s, and real jor real p and s. 

(ii) W is analytic in the domain Re p > 0; Re s > 0. 

(iii) W s -W+. 

(iv) For any (p , s ) with Re p = Re s = 0, which is a regular 
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point of W, poles of W(p , s) and W(p, s ) are simple and restricted to 
the imaginary s and p axes respectively. 

(y) dW/dp and dW/ds are positive semidefinite Hermitian for Rep = 
Re s = 0, except at poles. 

The proof of this theorem can be found on p. 33 of Ref. 8. We will 
interpret the above conditions on a physical basis. Assuming that a 
network realization consisting of reactances in the p and s-planes 
exists for W(p, s), condition i is fairly obvious, since the general loop 
impedance will be a real rational function in p and s. Condition Hi is 
also an obvious consequence of this reason, since the substitution of 
—p and — s in W{p, s) is equivalent to changing the sign of all ele- 
ment values and hence of every branch and loop impedance. Under 
the assumption of existence of a two element kind of reactance net- 
work corresponding to the given W(p, s), condition iv is also clear, 
because p is fixed as a pure imaginary number, the p-type elements 
can be considered as "frequency insensitive reactances," and their 
presence in a network consisting of pure reactances in the s-plane can- 
not create poles off the imaginary s axis. Similar reasoning justifies 
condition v for s fixed at any imaginary number, the positive semi- 
definiteness of dW/dp can be considered as an extension of the posi- 
tive slope of a reactance function in the one variable theory. 

The necessary and sufficient conditions for a two variable reactance 
function are not discussed separately, since scalars can be considered 
as a special case of a reactance matrix. 

If W{p, s) has a pole at p = p 0) independent of the value of s, po 
is said to be an s-independent pole of W. The following theorem (see 
p. 34 of Ref. 8) concerning such poles is important for the synthesis 
method to be given. 

Theorem 2: A two variable reactance matrix Wo(p, s) can be de- 
composed as 

W (p, s) = Wi(p) + W 2 (s) + W(p, s) 

where Wi and W 2 are reactance matrices in p and s, respectively, and W 
is a two variable reactance matrix with no p-independent or s-independent 
poles. 

Any given two variable reactance matrix W (p, s) , by virtue of the 
above theorem, can be realized as a series connection of networks hav- 
ing Wi, W2, and W as their impedance matrices, as shown in Fig. 1. 
Since Wi and W 2 can be realized by existing techniques (See chapter 
7 of Ref. 11) the given W can be realized if a method of synthesis is 
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Fig. 1 — Interpretation of theorem 2. 

found for W. Henceforth we assume that the given reactance matrix 
has no p-independcnt or s-independent poles. 

III. SYNTHESIS OF TWO VARIABLE REACTANCE MATRICES 

Let us assume that there is a passive n-port network representation, 
consisting of p- and s-type reactances, gyrators, and ideal transform- 
ers, for a given two variable n X n reactance matrix W{p, s). In such 
a network it is always possible to replace each s-type capacitor by a 
gyrator-s-type inductor combination and then isolate all the s-type 
inductors, of which we assume there are k, as shown in Fig. 2, without 
changing the impedance seen at the prescribed ports. If we further 
assume that the (n + fc)-port coupling network, consisting of p-type 
reactances, ideal transformers, and gyrators has a Z matrix, then the 
impedance matrix W(p, s) seen at the first n ports is given by 

W(p, s) = z u (p) - z i2 (p)[z 22 (p) + s/,r'z 21 (;;) (1) 

where Z(p) , the impedance matrix of the coupling network is given by 



Z(p) = 



Zu(p) Zizip) 

L*2i(p) Z 22 (P)J 

Since the coupling network is a lossless network in the p-plane 

Z = -Z f 
and we have 

W(p, s) = 2„(p) + z 12 (p)[z 22 (p) + sl k )~ l zUp)- (4) 

Next we show, by algebraic means, that every two variable reactance 

matrix can be decomposed into the form in equation (4), such that 



(2) 



(3) 
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Fig. 2 — Extraction of S-type inductors. 

Z(p) of equation (2) describes a lossless network. Once such a decom- 
position is found, we can realize the given W(p, s) by realizing Z(p) 
by any of the existing techniques (see chapter 7 of Ref. 11) and 
terminating it at its last k ports with unit inductors in the s-plane. 

To establish that any given two variable reactance matrix W(p, s) 
can be decomposed as shown in equation (4), we first expand W(p, s) 
and the expression on the right side of equation (4) about s = oo and 
find the expressions that relate 2 n , z-i 2 , and z 22 with the expansion 
coefficients of W{p, s). We then show that a set z llt Zi 2 , and z 22 , which 
satisfies the above relations and at the same time guarantees that the 
Z(p) of equation (2) is a reactance matrix in p, can always be found. 

The given two variable reactance matrix W(p, s) can be assumed 
to have no p-independent or s-independent poles by virtue of Theorem 
2 and hence can be written in the form 



+ B r (p) 
+ a T (p) 

where the Bi(p) are real polynomial matrices in p and the scalar 



= B n ( P )s + B&y- 1 + 

{p > S) a (p)s + afoW 1 + 



(5) 



g(p, s) = a (p)s r + a 1 (p)s r 1 + • • • + a r (p) (6) 

is the least common denominator of the entries in W(p, s). For any 
ordinary value of p, W(p, s) can be expanded in the neighborhood of 

s = oo as 9 



w(p,8) = A- l ( P ) + z4#- 

1=0 s 



(7) 



Expanding the right side of equation (4) in the neighborhood of s = oo 



i t 

ZtoZov.Xt 



*»(p) + *.(p)fe»(p) + sl^zUp) = Zll + (-1)' £ 51^512. ( 8 ) 

1=0 S 



For the equality in equation (4) to hold, we identify 
z u (p) = A_,(p) = W(p, oo) 



(9) 
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and 

A,(p) = (-1)Wm*I. I = 0, 1, 2, • • • . (10) 

Since the Z (p) formed out of *u , z 12 and z 2 2 

z(p) = r 2n t 2i2 i en) 

_ Z12 #22- 

has to describe a lossless network in the p plane, we must have 

Z = -Z f 
as given by equation (3) , and hence 

*u = -z» (12) 

and 

z 22 = — Z22 • (13) 

With the identification in equation (9), equation (12) is always satis- 
fied, since by equation (9) 

211 = W(p, «) = -W(-p, -co), 
and thus 2 U is uniquely determined. The problem is to chose a pair 
Z12, 222 to satisfy equation (10) and at the same time guarantee that 
equation (11) describes a lossless network in the p-plane. For Z(p) to 
describe a lossless network, it must be positive real and satisfy equa- 
tion (3). 

Before proceeding further, we would like to know more about A t {p) , 
the expansion coefficients in equation (7) . By equating the right sides 
of equations (5) and (7), 

B (p)s + B.W 1 + • • • + B r (p) 
= [ooCpK + ihtpW- 1 + • • • + a r (p)][A_ 1 (p) + g ^#J- (14) 

Equating coefficients of like powers of s on both sides of equation (14) , 
(seep. 207 of Ref. 12, Vol. II). 

aoWA-te) = B (p) 

Ox^A-xCp) + a (p)A (p) = 5 x (p) 

a 2 (p)A_i(p) + ax(p)A (p) + a (p)A 1 (p) = B 2 (p) ^ 
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ar(p)A-t(p) + Or-i(p)A (p) + ••• + a (p)A r -i(p) = B r (p) 

and 

a (p)A i (p) + a 1 (p)A i - 1 (p) + • • • + a r (p)A,._ r (p) = n for i ^ r. 

From equation (15) an expression for A t (p) can be written* in the con- 
venient form (see p. 14 Ref. 9) 



(16) 







B (p) 


Oo(p) 





■ 






B x (p) 


fll(p) 


a„(p) 
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(-1)' +1 

" «o +2 (P) 
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a 2 (p) 


a,(p) a„(p) 









B,(p) 
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a t+1 (p) 


a,(p) a,_,(p) • 


• fll(p) 






I 


= -1, 


0,1,2, ... 








B t 


= o n 


for I > r 








a, 


= 


for I > r 





where the (Z + 2) X (Z + 2) determinant is expanded formally in 
terms of its first column. In equation (16) the B\ are matrices, the a { 
are scalars, and the determinant is not a determinant in the usual 
sense. From equation (16), it can be seen that A x {p) is of the form 



A l {p) = 



real polynomial matrix in p 
a ' +2 (p) 



(17) 



Another important property of the Ai(p)'s is obtained from the rela- 
tion 

W(p, s) = -W(-p, -a) 
which implies 

a-ap) + i:4# = -a-x(-p) - z (-d' +i -4^- as) 

J=0 o * 

Hence by equating like powers of s 



* Alternate methods of obtaining these Ai(p)'s are by diffe rentiation of W(p, s) 

d l+1 W(p, s) 



Ai(p) = 



d8 l 



or by long division. 
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A l = (-l) l A\. (19) 

If, for the purpose of choosing a pair Z12, 222 that satisfies equation 
(10) and, at the same time, guarantees that the Z(p) of equation (11) 
describes a lossless network in the p-plane, we define Pdp) as 



r 



2i 



Ptlp) = 



Then 



PPi = 



Z 12^22 
^12^22 

I 

_Z- i2 Z 22 _ 



t t t ., 2t , 2t 

Z\2Z\2 212Z22Z12 2x2Z22Zi2 

f f * >7 2t ^ t 

2l2Z22Zl2 Z12Z22Z22Z12 2 12 2 2 2Z22 Z 12 

2 t 2 t t 2 at t 

2i2222Z]2 212222Z22Z12 Z\lZ 2 %L 22 L 22 



(20) 



it t ■ 

• • 212Z22Z12 

it t 

• • Z12Z22Z22Z12 

2 it t 

• • 212222Z22Z12 



It I t t I 2t t I It t 

LZl2222Zl2 Z12Z22Z22Z12 2]2222Z22 Z 22 ' ' " Z 12 Z 2 2Z21 12 - 



(21) 



In the above matrix, the entry in the zth row and ;th column is Zi2Z2 2 z 22 z,., 
and by equation (13) 



Z12Z22Z22Z12 ~ \ -U Z12Z22 Zj2 • 

Since we wish the equality in equation (10) to hold 

z l2 zi 2 ^iW\2 = (-i)'zi 2 z 2 2'zi2 = (—iyA i+i 

If we define Ti{p) as 

A (p) A t (p) A M (p) 

-A v {p) -A 2 (p) -A 3 (p) 

Ti(p) = Uv) Mp) AM 



(22) 
(23) 



-A l+l (p) 
A l+2 (p) 



_(-l) l A t (p) (-l) l A l+l (p) (-l)% +2 (p) ••• (-l)%i(p)_ 

(24) 
from equation (23) , we can see that 

T l = P,P; . (25) 
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Equation (25) suggests that a way of obtaining a pair z 12 , z 2 i would 
be to form the matrix Ti(p), factor it in the form of equation (25), 
and then try to identify z ]2 and z 2 2 from these factors. We do not know 
in advance if the matrix Ti(p) formed from the expansion coefficients 
of W about s = oo can always be factored as indicated in equation 
(25) ; hence we first study the properties of Ti(p), to see if it can be 
factored in the desired form. 

Consider the matrix Ti (p) when I = r,r being the s-degree of g (p, s) , 
as given in equation (6) , 



T r = 



A 


Ai 


A 2 


4,-. 


A r 


-4, 


-A, 


-A, 


-A, 


— A r+ i 


A s 


A a 


A, 


A r+X 


A r +2 



(-D-Mr-, (-D-^r (-l)-'A r+1 ••• (-l)^A 2r -l (-l)^A 2r -l 
|_(-U r i4r (.-D'Ar+l (-D r A r+ 2 ••• (-l) r A 2r -l (-lYA ir 



• (26) 



The matrix obtained by deleting the last column and row in equation 
(26) is T r -i, and by equation (15) it is easy to see that the last column 
is a linear combination of the first r columns. Hence* 

rank T T = rank T 7 ,., 

and 

rank T t = rank T r -, for I ^ r - 1. 

The rank of T r _i is connected with the s-degree 8 8 [W(p, s) ] of W(p, s) 
which is defined in Definition 3 (see p. 10 of Ref . 9) . 

Definition 3: The s-degree of a rational two variable matrix W(p, s) 
is obtained from the rule 

s = degree of W (p, s) = S,[W(p, s)] = max 8[W(p , s)] 

where 8[W(p , s) ] is the McMillan degree (see part II of Ref. 13) of 
W(p , s). For any fixed p , W(p , s) is a matrix of rational functions 
in s with its McMillan degree uniquely specified; hence the above 
definition uniquely specifies the s-degree of W(p, s). The relationship 
between the s-degree of W(p, s) and the rank of TV-i is stated formally 
in the following lemma. 

Lemma 1: The rank of r r -i(p) *s equal to the s-degree of W(p, s). 



*By the rank of rational or polynomial matrix we mean the "normal rank," 
which is defined to be the rank everywhere except at a finite number of values 
of the variable. 
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The proof of this lemma for the one variable case can be found in 
Ref. 14 and on p. 200 of Ref. 10, and for the two variable case on p. 17 
of Ref. 9. 

To show that the matrix Tr-i(p) can always be factored in the form 
of equation (25) , we need the following lemma. 

Lemma 2: The matrix T T - X (p) defined by equation (24) for I = r — 1 
satisfies 

(ii) T T ^(ju>) is Hermitian and positive semidefinite. 

Proof: 

Since A, = (-1)'A[ by equation (19), the proof of i is readily seen 
from equation (26) 

T r . x = Tl, . (27) 

To prove (ii), we first notice that by Theorem 1, for any real <o, W(ja>, s) 
has only simple poles, which are restricted to the imaginary axis in 
the s-plane. Hence W(ju, s) can be expressed in the partial fraction form 

ire*., it = a-CM + t ^gk) (28) 

where, Ri(j<a) are the residue matrices at the poles jai(u>), and the a,-(w) 
are real. 

It is shown in Appendix A that the fl.O'co) are Hermitian and positive 
semidefinite for each co. Now, if each term in the sum on the right side 
of equation (28) is expanded about s = °o , we have 

W(j», s) = A^fr) +EE *S£ «<(*)■ ( 29 ) 

For the purpose of comparison, equation (7), written with p = jw, is 

W(p,,s) = A- J (jo>)+ Z4#- (30) 

= * 

The right sides of equation (29) and (30) are expansions of W(ju, s) 
about s = c» , and because of the uniqueness of a power series expansion 

AM = ECjttO^Gw). (31) 

i = i 

By noting that the a, are real and the R t (ju) are Hermitian and positive 
semidefinite for each co, we have 
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Aoiju) = EB,G») ^ (32)* 

i = i 

A,(j») = j E«*fi«(iw) (33) 

»-i 

A,(jw) = - E c&ltfr) ^ (34) 



A 4 ,„_ 3 (; W ) = -j £ «?- 3 Ri(p>) 

i-1 

A im - 2 {jo>) = -J2 etr-Rtfr) ^ 
t=i 

Antics) = -i i: as- 1 ^) 

i=i 

4 4 *(jw) = E cfRtiju) ^ 0. 



(35a) 
(35b) 
(35 c) 
(35d) 



By direct substitution of equation (33) into equation (24), T r -i(jo)) 
can be written as 



-V-i(M = E 



fl, 


jotiRi 


-at'Ri 


• • G*)" 1 ^ 


-jafii 


a'Bi 


ja'Ri 


• • -(*,)'«, 


-uBi 


-jafii 


«Jfl. 


0a) r+ '«. 



(36) 



The matrix sum on the right side of equation (36) can be written 

~R { • • • 0" 



Tr-lW = £ U 



••• 
••• 

• • • 0. 



L* 



(37) 



where 



* By the notation A ^ or A ^ 0, we mean that the associated Hermitian form 
of A is positive semidefinite or negative semidefinite. 
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Li = 



l n • • • 

j«,ln In • • • 

a*l n 1„ ••• 

L(-i) r (i«) r_1 i„ o o ••• i. 



(38) 



Since each Ri(ju>) is Hermitian and positive semidefinite for each w, 
the sum on the right side of equation (37) is also Hermitian and positive 
semidefinite. Hence, we have proved the lemma. 

We have shown that Z\._ x , a matrix of rational functions, is para- 
Hermitian and positive semidefinite on the imaginary axis. Such a 
matrix can always be factored in the form shown in equation (25), 
(see p. 133 of Ref. 15). It is tempting to factor T r - X at this stage and 
find z 12 , z 22 to satisfy the required conditions, but we will factor a\ r T T - X 
instead of T r _i for the reason that the factors would be polynomial 
matrices. 

From equation (17) we can see that a 2 r T r _i is a polynomial matrix 
in p. To be able to factor a 2 r T r - X in the required fashion, we have to 
show that t = a 2 Q r T r - t is para Hermitian and positive semidefinite on 
the jco axis. To do this, we obtain the required additional information 
about the polynomial a (p) from the following theorem. Since the 
theorem contains more information than we need at this point, we will 
only state it here; a proof is given in Appendix B. 



Theorem 3: If 



w<p,.)- B <W + B >W 



+ B r (p) 



a () (p)s r + a^py + • • • + a r (p) 



is a two variable reactance matrix, then for all i = 0, 1, • • • , r 
(i) '—■ is a reactance matrix in p 



di 



(n) 
{Hi) 

(*) 



a, 



is a reactance function in p 



a, has all its zeros on the ju axis and these are simple 

XBiX for all constant real n X 1 vectors, X, is a reactance 
XB i+l X function in p 
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From this theorem, a„(p) can be represented as 

flo(p) = V II (V + «!) - ±o«(-p) (39) 

where y = or 1. Hence 

h-l = «oTr-X = t r f _, . (40) 

From the form of a shown in equation (39) and Lemma 1, it can be 
seen that 

*V-iG») ^ o (41) 

except when simultaneously, v = 1 and r is odd; in which case 

franco) ^ 0. (42) 

We will assume that f r -i(j<a) ^ in developing the synthesis procedure 
and discuss the needed modification when T r -i(ju>) ^ later. 

If the s-degree of W(p, s) is equal to k, by Lemma 1 the rank of T r - X (p) 
and hence of T r -M is k. Since f,^ = tl_ t and t r -i(ja) ^ there 
exists a factorization 18,17 

f r _M = M(p)W(p) (43) 

where M(p) is annrXfc polynomial matrix and has a left inverse M _1 (p) 
which is analytic in Re p > 0. 

From the definition of t r -i, we have 

do 
M{p) can be partitioned into n X k blocks M { (p) 

M (p) 



(44) 



M(p) = 



M,(p) 



MUph 



(45) 



and hence 



m\ p ) = \ml(p) ; mI(p) : • • • ; idfo)]. 



(46) 



Now by comparison of equation (45) with equation (20), we can im- 
mediately identify a suitable z 12 as 



Z\o 



M Q 



(47) 
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To find a suitable z 22 , if we define T d as 
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-A, 


-A 2 


-A 3 ■ 


■ -A T 


A 2 


A 3 


A 4 ■ 


A r+l 


-A s 


-A, 


-A 5 ■ 


■ -A r+2 



tm-= -~> ; 14 r . " (48) 

_(-l) r A r (-l) r /l r+1 (-l) r A r+2 ■■■ (-l) r A 2r J 
from equations (20), (21), and (25) we see that z 22 must satisfy 

~ Mz 22 M f = T d . (49) 

a 

Even though equation (49) does not uniquely specify z 22 , we can choose 
for z 22 

z 22 = al'M^TM' 1 *. (50) 

From equation (19) and the definition of T d , we see that T d = - T (/ 
and hence 

Z 22 = ^22- » ' 

We now notice that by construction, the pair z 12 , z 22 defined by 
equations (47) and (50) satisfies 

(10) 



{-l) l z l2 z l 2 ,z\ 2 = A , 

for all ^ I S 2r - 2. Our aim is to find z„ and z 22 that satisfy equation 
(10) for all 1^0. It is not immediately clear that the pair z 12 , z 22 defined 
by equations (47) and (50) satisfy equation (10) for all U 0. 

To see that the chosen pair z v2 , z 22 does indeed satisfy equation (10) 
for all I ^ and not just for g I g 2r - 2, we introduce the generalized 
companion matrix fi(p) defined by 1 



b(p) = 



n 

o n 



.&1 _ 



In 

o„ 

flr-1 



o„ 

1„ 



«0 



0, 



_0, ln _«!! 

a a J 



(52) 
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From equation (15) it can be seen that 

T d = -T T -,a. (53) 

Hence, by equation (43) 

z 22 = -aTM^Tr^QW^ 
= -Af _1 JlfM t OM" lt 



t^Tur-lt 



and by equation (51) 
Hence 



= -MQM 



z 22 = -MW-'* = M _I fi + M. (54) 



and 



2 2 22 = -AT'O^/M^M"" 1 ' 
= -M -1 X2 2t M 



4 = (-l)' -1 iirV ilf; Z > 0. (55) 



From the definition of Q, we see that </(p, £) is its minimal polynomial, 
and hence the matrix polynomial 

g(p, Q) = a Q r + a,Q r_1 + • • • + a r l n r = nr (56) 

and hence 

g(-p, fit) = nr . (57) 

By equation (55) 

g(p, z 22 ) = il/- 1 (-l)'" 1 ^ + i-iy-W-^ + • • • + a,l nr M, 
and by equation (57) and Theorem 3 

g(p, z 22 ) = ± ar'foC-p, tf)]M = 0*. (58) 

From the last equation in equation (15), from equation (58), and from 
equation (10), which holds forV^ I ^ 2r - 2, 

= — 3, 2 [ai222 -2 + a 2 222~ 2 + • • • + flwSi^Zw 

— __ 2r_1 f 

^0^12^22 Z]2 . 
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Hence 



and by induction 



A 2r _, = (-l)' r " 1 «i«a£r 1 z t 



12 > 



71/ — (— 1) 212222Z12 



for all Z ^ 0, which is the same as equation (10). 

We thus have a set of three matrices * U) z 12 , z 22 such that the in- 
finite set of equations obtained by equating the right sides of equations 
(7) and (8) are satisfied. Hence the right side of equation (4) and W(p, s) 
have the same Taylor's series expansion in the neighborhood of s = 00 . 
By analytic continuation, for all p and s 

W(p, S) = Z n (p) + *ia(p)fen(p) + Sl,] _1 z! 2 (p), 

where z u , z 12 , and 2 22 are defined by equations (9), (47), and (50), 
respectively. 

We have thus succeeded in decomposing W(p, s) as shown in equation 
(4). It now remains to show that Z(p) formed from the chosen z u , 
2 12 , and z 2 2 

I" 2n(p) 2ia(p)l 

_-Z f n (p) 2 22 (p)_ 



Z(p) = 



(59) 



TF(p, 00) 

Ml(p) 
L a^Cp) 



MM 

«o(p) 



ar(p)M-\p)T d (p)M- l \p) 



is a reactance matrix. 

To show that the Z(p) in equation (59) is a reactance matrix, we may 
choose any standard test, but we will choose the one given below since 
it is particularly suited for the problem at hand (see pp. 117 and 123 
of Ref. 11): 

Lemma 3: The necessary and sufficient conditions for a square matrix 
Z(p) to be a reactance matrix are: 

(i) Z is rational and real for real p. 

(ii) Poles of Z(p) are simple and restricted to the imaginary axis. 
(Hi) Z + Z f = 0. 
(iv) Residue matrices are positive semidefinite- Hermitian. 

Since all the entries of Z(p) hi equation (59) are real and rational, con- 
dition i is satisfied. 
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From equations (13a) and (15), z u = A- r = B /a is a reactance 
matrix by Theorem 3; hence its poles are simple and restricted to the 
imaginary axis. Also, the pole at p = <x>, if any exists, is simple for 
this block. Since M is a polynomial matrix, it is clear that the poles of 
the off diagonal blocks z 12 and — z J2 are in the zeros of a and hence by 
Theorem 3 the poles of z 12 are restricted to the ju axis. However, it is 
not clear that these poles are simple. To show that they are indeed 
simple we will use the fact that A denned by 

A = a 2 A (60) 

is a polynomial matrix. From equations (10) and (47) 

A. = «g (61) 



oS" 1 J • 



(62) 



and from equation (39) a = ±a . We first consider a = al in which 
case 

-^0 — r-1 

_flo J I 

Equation (54) then shows that A = A<J and A (ju>) ^ 0. Hence there 
exists annX? polynomial matrix, Q, such that 

A = QQ+ (63) 

where q is the rank of A . Equations (62) and (63) are two different 
factorizations of A , hence: 17 

^ = Q[K\0 kx(T - q) ]V (64) 

where V(p) is a k X k para unitary matrix, that is, FV t = l k . Since Q 
is a polynomial matrix, and V(p) being para unitary can have no poles 
on the imaginary axis (see p. 186 of Ref. 11), the left side of equation 
(64) can have no poles on the imaginary axis. Hence a T Q ~ l , which has 
all its zeros on the jta axis, must divide M . Thus z 12 has all its finite 
poles in the zeros of a . By Theorem 3, the zeros of a are simple and 
restricted to the jta axis. In the above, we have assumed that a = a ; 
if a = — a and r is odd, the same proof holds; if r is even we can con- 
struct a similar proof by factoring — A instead of A . 

To show that the pole of z 12 at p = <x>, if any, is simple. Consider 
the following representation for A obtained from equations (15) and (17) 

. _ a Q B 1 — c^Bq _ B± <h _ B Oi . . 

al a v a a a 
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Since B x /a x and B /a are reactance matrices and a x /a is a reactance 
function, according to Theorem 3, the right side of equation (65) be- 
haves as Kp" near p = °o , where K is a constant matrix and v is an 
integer such that — 2 ^ v ^ 2. But 2 12 satisfies 

and hence the pole of z 12 at p = » , if any, must be simple. 

We now have to show that the z 22 block also satisfies condition (it) 
of the lemma. By equation (54) 

z 22 = -WaW- 1 = M-WM. (54) 

Since M~ l is analytic in the open right-half plane and £2 has all its poles 
in the zeros of a , by equation (54) the poles of z 22 are restricted to the 
ja axis. To show that these poles are simple we will prove by contradic- 
tion that a z 2 2 is polynomial. 

From equation (52), the definition of 12, and equation (54) we see that 
if a Q z 22 has a pole of multiplicity a at p = ja> a . In the neighborhood 
of this pole, we have the approximation 

K 



(P - J«o) 
where if is a constant matrix and a is a positive integer, and 



(66) 



alzl, » ; 7—7: ■ ( 67a ) 



(p - i^o) 2 " 

Now by equation (55) z 22 = — M -1 fl 2t M , and hence in the neighborhood 
of p = j«o 



n - 2 - ' M (67b) 



a z 



(p - &o) 



where K l is a constant matrix and /3 is a positive integer. Since the poles 
of a z 22 are contained in the poles of M' 1 , ^ 2a. By comparison of 
equations (67a) and (67b), which must be equal, it is clear that either 
a = p = or Ki = K 2 = 0. Since z 22 = —z 22 , K = K*, and hence 
K 2 = KK* = implies that K = 0. Thus a 2 2 2 can have no poles on 
the ju axis and this, coupled with the fact that z 22 can have poles only 
on the jco axis, guarantees that a z 22 is always polynomial. We therefore 
conclude that all the finite poles of z 22 are in the zeros of a , and their 
multiplicity cannot exceed that of the corresponding zeros of a . Hence, 
again by Theorem 3, all the finite poles of z 22 are simple and restricted 
to the jo) axis. 
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To show that the pole at p = oo of z 22 , if any, is simple, consider 
equation (15) written in this form: 



(68) 



a 








4 _ ^i.*i _ ^o.^i 
a, a a «o 








a 2 aj a «o 


fe.^i. 


_ Bo (h 


CL 2 Cb\ -Do 


_«i Oo 


a a _ 


di do a 



Owing to the reactance nature of B,/a, and ai/a i + l by Theorem 3, 
and from the form of A { shown in equation (68), near p = <x> , A t behaves 
as 

A { & K iV Vi (69) 

where K { is a constant matrix and v { is an integer such that 

i+ 2 ^ v, ^ -(i + 2). (70) 

Also from equation (10) 

*u4r*Li = (~1)% ~ ±K iV vi . (71) 

Since 2 J2 has at most a simple pole at p = oo , in the neighborhood of 

p = go 

«,, « tfp' (72) 

where K is a constant matrix and I is an integer such that I ^ 1. If 
z 22 behaves as K 22 p m near p = oo , where .K22 is a constant matrix and m 
an integer, then by equation (70), (71), and (72), (i + 2) ^ im -+- 21 ^ 
— (i + 2) . For such to be true for any fixed I and all integral i ^ 0, 
m has to be less than or equal to unity. Hence the pole of z 22 at p = 00 , 
if any, is simple. 

We have thus shown that condition ii of Lemma 1 is satisfied for 
each block in Z(p), and hence Z(p) also satisfies it. 

Since z xx is a reactance matrix, z u = — z n and2 2 2 = — z 22 by equation 
(51), we have Z = — Z f and thus condition (in) of the lemma is also 
satisfied. 

Now to complete the proof that Z(p) is a reactance matrix, we have 
to show that the residue matrices at the poles are positive semidefinite 
Hermitian. To do this we need Lemma 4, which follows from the defini- 
tions of a two variable positive real and two variable reactance matrices 
(see p. 34 of Ref. 8). 
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Lemma 4. Ij W(p, s) is a two variable reactance matrix with no p-in- 
dependent or ^-independent poles, W[p, s(p)] is a reactance matrix in p 
jor any reactance junction s(p). 

To prove that Z(p) satisfies condition iv of Lemma 1, which requires 
that the residue matrix of Z(p) at any of its simple poles on the jcj axis 
is positive semidefinite Hermitian, we note that at any pole, p = ju>, 
of Z(p), if we set 

2lp 
s(p) = 2 , 2 for « < 00 

p + Oio 

= Ip for co = °o 

in 

W(p, s) = 2,, + 212(222 + s^) _1 z, 2 

[which is equation (4)] then by Lemma 4, W[p, s(p)] is a reactance 
matrix in p for all positive I. Since Z(p) is real for real p and Z = — Z + , 
the residue matrix H at the pole p = ju is Hermitian; if we write it as 

H = //n Hyl (73) 

_xi 12 ti 2 2_ 

then, K, the residue matrix of W[p, s(p)] at p = jo> is given by 

K = H u - //,,(// 22 + ZUr'H* . (74) 

Since H, H u , and H 22 are Hermitian, there exist unitary matrices f/ x 
and U 2 such that 

A n = U*//nt/, = diag [rf, , d 2 , ■ ■ ■ , d n ] (75) 



and 



Hence 



where 



A 22 = V*H 22 U 2 = diag [X, , X 2 , • • • , X,]. (76) 

\J*KU 1 = A,, - J 12 (A 22 + Ih)- 1 ]* (77) 



J 12 = \J*H i2 U 2 . (78) 

If J 12i denotes the ith column of J 12 , the right side of equation (77) can 



be rewritten as 



JJ*KU, = A n - J2 ^—^ J 12 . Jf 2 . • (79) 
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Since K is the residue matrix of a reactance matrix, for all I > 0, K is 
positive semidefinite. A„ is also positive semidefinite, since H lt is 
the residue matrix of the reactance matrix z u . J 12 .Jf 2 ,- is obviously 
positive semidefinite and the left side of (79) can be positive semi- 
definite for all positive I only if all the X,- are nonnegative. Hence A 22 
and #22 are positive semidefinite. 

To show that H is positive semidefinite, we will show that H' defined by 



iv = (u* + u?)ff(tfi 4- u 2 ) = 



(80) 



An Jn 

_Jl2 A 22 J 

is positive semidefinite. For this purpose, consider the Hermitian form 



[XI X* 2 ] 



A„ 

Jf 2 



D, 






= X*A u Xi + X* 2 Jf 2 Xi + X*/ I2 X 2 + X* 2 D 22 X 2 



(81) 



where 



Since 



D 22 = A 22 + ili 



I > 0. 



JJ*KU 1 = A u - J 12 D-IJ* 

is positive semidefinite, we obtain from equation (81) the following 
inequality: 



[XI XI 



A„ 

LJ?2 



■A: 

Do. 



X 

x 2J 



^ X* 1 J i2 D- 2 1 2jf 2 X 1 + X\J 12 X 2 + X* 2 Jf 2 3Ti + X* 2 A 22 X 2 



(82) 



Since the right side of equation (82) can be expressed as G*G, where 
G = [D l 2 / 2 2 Jf 2 X 1 + An'XJ, the Hermitian form in equation (81) is 
positive semidefinite for all I > 0; by a continuity argument we can 
see that H' and consequently H are positive semidefinite. 

We have thus shown that Z(p) does indeed describe a lossless network 
in the p-plane and thus W(p, s) has the network representation shown 

in Fig. 2. 

In the development of the synthesis procedure we assumed that 
ao r O'a))T r _i(ia)) = tr-i(jv) ^ 0. If simultaneously, a (p) is an odd 
function of p [in other words v = 1 in equation (39)] and r, the s-degree 
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of the least common denominator of W(p, s) is odd, then T r -i(ja>) ^ 0. 
In this case we factor — t r -i(p) which is para Hermitian and positive 
semidefinite on the jta axis. We will then have 

-TV-, = MW (83) 

and hence, as before equation (44), 

T -MM! 

r_1 ~~ n r '* 

It is then clear that the identification of z l2 and z 22 can be done in 
exactly the same way as when T r - X (ju>) ^ 0. 

It is of importance to notice that the number of s-plane inductors 
used in the realization of Fig. 2 is equal to the s-degree, 8 8 \W(p, s)] 
which in general is smaller than the number required in Koga's tech- 
nique. Appendix C shows that S s [TF(p, s)] is the minimum number of 
s-plane inductors required in any realization, and that if a realization 
is minimal in the variable s it is automatically minimal in the variable 
p, the minimum number of p-type reactances needed in any realiza- 
tion being the p-degree, 8 p [W(p, s)]. 1B 

The main result of this section can be conveniently put in the form 
of a theorem : 

Theorem 4: Every two variable reactance matrix W(p, s) can be realized 
as the impedance seen at the first n-ports of a lossless (n + k)-port con- 
sisting of 8 p [W(p, s)] reactances in the p-plane, terminated at its last k 
ports with 8,[W(p, s)] unit inductors in the s-plane. Furthermore, such a 
realization uses the minimum possible number of reactances of each kind. 
(The roles of p and s are completely interchangeable.) 

Since several of the proofs involved in establishing Theorem 4 were 
rather indirect and lengthy, while the procedure for synthesis, sum- 
marized in Section IV, is itself rather simple. 



IV. SUMMARY OF SYNTHESIS PROCEDURE 

Given an (n X n) two variable reactance matrix W (p, s) , decom- 
pose it as 

W (p, s) = W l (p) + W 2 (s) + W(p, s) 

where Wi and W 2 are reactance matrices in p and s, and W(p, s) is a 
two variable reactance matrix with no p-independent or s-independent 
poles. Such a decomposition is always possible by Theorem 2. 
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Expand W(p, s) as 

W<p,8)mA- l ( P )+ Z4# 

J-0 s 

[which is the same as equation (7)] where the A(p)'s may be ob- 
tained by equations (16) or (16a) or by long division. 

Find g(p, s), the least common denominator of the entries in 
W(p, s) and express it in the form 

g(p, s) = a (p)s r + a,(p)s r_1 + • • • + a {p). 

[which is the same as equation (6) ] . 
Form the (nr X nr) matrix 7V-i (p) , defined by 

A (p) A x (p) A 2 {p) ... A r _,(p) 

-A,(p) -A 2 (p) -A 3 (p) ... -A r (p) 

Mp) A 3 (p) A<(p) ••• A r+1 (p) 

U-ir'iUifo) i-i)" x A r (p) (-iy- l A r+l (p) • •• (-i) , " 1 A, r . 1 ^)J 

which is equation (24). 

Factor 7\._i(p) = al T T r _ x (p), a polynomial matrix, as 

t r -M = MW (43) 

unless simultaneously, a = — a and r in equation (6) is odd, in which 
case factor — T r -i(p). The factorization must be such that M is a 
(k X nr) polynomial matrix with k = rank of T r - t (p) and M~\ the 
left inverse of M analytic in the open right plane. The existence of such 
a factorization is guaranteed by Lemmas 1 and 2. 
Partition M(p) into (n X k) blocks of equation (45) 

M (p) 

M,(p) 
M(p) = 

Form the (nr X nr) matrix fi(p) defined by equation (52) 



REACTANCE MATRICES 187 

With the identification of equations (9), (45), and (54)* 
z u (p) = ,4_,(p), 

, M (p) 

and z x2 = r~r~ > 
a (p) 

and z 22 = M~\p)^(p)M(p), 
the decomposition 

W(p, s) = Zll (p) + z l2 {p)[z 22 {p) + sl t ] _1 z; 2 (4) 

is obtained. Notice that this is equation (4) . It should also be noticed 
that W{p, s) can be decomposed as in equation (4) even if it has 
s-independent poles, since the assumption that W(p, oo) is finite is 
enough to guarantee the validity of the procedure. For network reali- 
zation it is usually more convenient to remove both p-independent 
and s-independent poles; we therefore removed them at the start of 
the procedure. 

To realize W(p, s) as the impedance of a passive network, we per- 
form the following operations. 

Form the (n+k X n+k) impedance matrix Z{p) of the coupling 
network 



Zip) = 



, , , MM 

A - i{v) Wf 

_M^) M - 1(p) Qt (p)M(p) 

- flo(p) 



(84) t 



Realize Z(p) as a lossless (n+k) port network in the p-plane and 
terminate its last fc-ports with unit inductors in the s-plane. Also 
realize the reactance matrices TFj (p) and W 2 (s) as lossless p-plane 
and s-plane n-ports, and connect all three networks in series as shown 
in Fig. 1. The given W (p, s) is thus realized as a passive network. 

V. AN EXAMPLE 

It is desired to synthesize the two variable reactance matrix:}: 



* Equation (54) is used to determine z-iip), in preference to equation (50) since 
equation (54) is easier to compute. 

t Equation (84) is the same as equation (59) except that for the 0m block 
equation (54) is used instead of equation (50) for the reason mentioned in the 
previous note. 

$ This example was given by Koga, (see p. 50 of Ref . 8) . 
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W Q (p, s) = 



(p 2 + W + 1) ps- 1 
(p + s)(ps +1) p 4- s 



ps — 1 ps + 1 

p + s p + s J 

Since W {p, s) has no p-independent or s-independent poles the first 
step 1 of Section IV need not be performed, and W (p, s) = W(p, s). 
The least common denominator of the elements of W(p, s) is 

g(p, s) = ps 2 + (p 2 + l)s + v 
[which is equation (6) ] , and hence 

«o(p) = P, ai(p) = (p 2 + 1), a 2 (p) = p, and r = 2. 

The least common denominator of the minors of W(p, s) is also g(p, s) 
and hence 

k = S.[W(p, >)] = 2. 

In the expansion, equation (7), 

w(p,8)-a- 1 <p) + z4# 

1-0 o 

by the formula of equation (16) or by long division 



A-i(p) = 



V 



p 2 + 1 p 
L V 2 P'-i 



Mv) = -3 



A l(P) = p3 



A 2 (p) = 



(p+D 2 pV + D 
Lp 2 (p 2 + i) p 2 (p 2 -i)J 
(p 2 + i) 3 pV + i)1 
LpV + i) pV - i)j ' 

p 8 + 3p 6 + 4p 4 + 3p 2 + 1 p'Xp 2 + 1) 

pV + i) pV - DJ 

r r _i(p) = 7\(p) defined by equation (24) is 

_p*(p»_J-l)S _p4( p 2_|_ 1 ) | p(p 2 + l)» P 6 (P*+1) 

-p 4 (p 2 +l) -p 4 (p 2 -l) ' p 6 (P 2 +l) p 5 (p 2 -l) 

_p( p 2_(_l)3 _ p 5(p24-l) I p8 + 3p6_|_4 p 4 + 3p2 + 1 p 6( p 2 + 1 ) 

I 
-P B (P* + 1) -P 6 (P 2 ~1) ' p 6 (p 2 + l) p 6 (p 2 -l) 



Tr-M = - 4 
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The polynomial matrix T t (p) = ptT^p) is factored by the method in 

Ref. 16 as equation (43) 

p(p 2 +l) -p(p 2 +1) 

PKp-i) -pKp+i) 

p 4_p3 +2 p 2 -p + l -(p<+p 3 +2p 2 +p + l) 



T x {p) = M(p)MKp) = (2)j 



p4_p3 



X 



-p(p*+l) -p 2 (p+l) 



-(p4 + p3) 

p*+p*+2p i +p + l p*+p l 



(2)* _ p (p2 + i) p'(p-l) -(p 4 -p 3 +2p 2 -p+l) ~(P 4 -P 3 )_ 

The (4 X 2) matrix M(p) is partitioned as equation (45) 
MM 



M(p) = 



-Milp). 



pV + i) -pV + i) 

vXv- i) .Jl^i?. 

p V - V + 2p 8 - p + 1 ~(p 4 + p 3 + 2p 2 + p + 1) 

v'-v -iv+v 3 ) 

To find M- 1 (p) , a left inverse of M (p) , it is enough to find a left in- 
verse of Mq if it exists, since 

[Mo 1 ! 0] 



Mo 
LMJ 



= 1* . 



In our example k = 2 and M is a nonsingular matrix and hence 
M -1 (p) is given by 



M _1 (P) = CT. 



-1 



-p 2 (p + i) p(p 2 + i) ! o o 



2p 3 (p 2 + 1) __ p2(p _ 1} p(p2 + 1} i _ 
From the definition of fi, equation (52) 



Q(p) = 









1 














1 


-1 





p 2 + l 

V 








-1 





p 2 + 1 
p J 
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Using equation (54) 



-lnh 



z 22 = M~ l QM 



p 2 + 1 
2p 



(P + I) 2 
2p 



(P ~ I) 2 ?> 2 + l 
2p 2p 

Hence the coupling network formed by p-type elements has the 
(4 X 4) matrix of equation (84) 



Zip) = 



A_, 


Mo 
a 






-*f M 

a 


" l Q f M 






p 2 + l 
2 P 


V 


p 2 + l 
(2)V 


p 2 + l 
(2)*p 


V 


V 


p- 1 

(2)* 


p + 1 

(2)* 


p 2 + l 


v + 


1 p 2 + l 


(P + I) 2 


(2)*p 


(2)* 


2p 


2p 


p 2 + 1 


v - 


1 (p - D 2 


P 2 + 1 


(2)*p 


(2) 


* 2p 


2p . 



Z(p) can be verified to be lossless, and the given W (p, s) can of course 
be realized as the impedance seen at the first two ports of Z(p) when 
it is terminated at its last two ports by unit s-plane inductors. 



VI. CONCLUSIONS 

The synthesis method for two-variable reactance matrices developed 
here, in general yields a nonreciprocal coupling network even when 
the given two-variable reactance matrix is symmetric, and if a 
reciprocal coupling network is desired, Koga's method for generating 
a reciprocal network from the nonreciprocal one can be used. 8 This 
procedure generally yields a reciprocal network at the cost of in- 
creased numbers of elements of both kinds. 

This method of synthesis of two-variable reactance matrices has 
been successfully applied to the synthesis of lumped-distributed RC 
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networks which are important in microelectronics circuits. 7 In prac- 
tice, the only laborious step in the synthesis procedure is the factoriza- 
tion of polynomial matrix in the desired form. Of great importance is 
the approximation of desired characteristics by rational functions in 
two-variables; any work in this area would greatly enhance the use- 
fulness of the two-variable theoiy. The synthesis problem of n-vari- 
able positive real functions, for which many applications can be 
found, 7 can be reduced to the synthesis of (n+1) -variable reactance 
matrices. 21, 22 when n = 1 the two-variable method developed here 
gives rise to a new method of passive RLC synthesis, which is no more 
complex than the existing methods. 
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APPENDIX A 

Partial Fraction Expansion of W(ju, s) 

To guarantee the factorization of T r -i(ja) as MM 1 " we needed Lemma 
2, which asserts that T r - X is para Hermitian and that T r .i(ju) ^ 0. 
In the proof of Lemma 2 we used the fact that 22 ,■(;&>), the residue 
matrices of W(ju, s), are positive semidefinite. The proof is given below. 

Under the assumption that W(p, s) has no p-independent of s-in- 
dependent poles, for each real co the s-plane poles of W(j<a, s) are simple 
and restricted to the imaginary s-axis by Theorem 1. Hence, for any 
fixed co, we can write W(jw, s) as 

W(ju , s) = A.,^)+t^f^ (85) 

where the a.-(co) are real and the 22, (co) are the residue matrices at the 
poles jciiiw). As in equation (9), r is the s-degree of g(p, s), the least 
common denominator of the elements of W. 

By complex conjugation on both sides of equation (85) 

Since W and A_ x are rational in jta, 

W*(ja, s) = W(-ja,, a*) 
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and 

Atiiju) = A-ii-ju). 
Hence equation (86) becomes 



W(-p>, s*) = A.^-jco) + £ 



R%(o>) 



fz[ S* + jotii 



and 



-WC-jco, -s*) = -A-x(-jw) + £ 



Rl(co) 



fri s* — #*,-(«) 
Since equation (88) is an identity for s*, we have 



-W(-jw, -s) = -A-^-ico) + £ 



R*.(co) 



(87) 



(88) 



(89) 



, =1 s — jotifa) 
From the definition of a two variable reactance matrix, 

W(j<a, s) = -W(-jo>, -s) 
and by equation (19) 

A_,0'ai) = -A-^-jco). 

Hence, by comparison of equations (85) and (89) we have the desired 

result 



^..(co) = R*(co). 



(90) 



To show that the Ri(o)) are positive semidefinite for each oj, we first 
notice that if 

W<p, s)=^\ 
g(p,s) 

where if/(p, s) is a polynomial matrix and g(p, s) is the least common 
denominator of the entries in W, 12<(w) in equation (85) is given by 
(see p. 308 of Ref. 19) 



7? /„\ _ M>>«) 

Ki{<j)) = —, r 

5s 



(91) 



P = I 01 



Denoting dg/ds by g, and dr^/ds by $, , for any ?i X 1 constant matrix 
X, (p. 39 of Ref. 8) 



X*R { X = 
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(X* *X)g. - g(X^,X) "|- 1 
(X*^) 2 _L, U 



[ 

L^ \x*Jx) J 



i-i«((»i 



p=|U 



Hence, if X"WX ^ 

X*RMX = [|- (X*T7X)- 1_ 1 ' . (92) 

\_OS Jp-l'w 

•-iai(u) 

From definitions 1 and 2, and Theorem 1, X*WX is a two variable 
positive function and for Re p = Re s = 

Re [X*WX] = 

and 

d 



(X*WX) ^ 0. 

uS 



Hence 



for 

Re p = Re s = 

and consequently the left side of equation (90) is nonnegative. 

Thus we have proved that the residue matrices, #<(w), are positive 
semidefinite Hermitian for each «. 

APPENDIX B 

Proof of Theorem 3 

Theorem 8: If 

w , * B n (v)s r + B l (py~ l + • • • + BM 
W(V > S) = a (p)s' + afrV 1 + • • • + a r (p) 

is a two variable reactance matrix, then for all i = 0, 1, • • • r 
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Bi/di is a reactance matrix in p. 
CLi/di+i is a reactance junction in p. 
a, has all its zeros on the j axis, and these are simple. 
XBiX/XB i+1 X for all constant n X 1 vectors, X, is a reactance function 
in p. 

Proof: For any constant n X I matrix X, 

X *WX = ( x * g o*)* r + (yg^r 1 + • • • + (X*B T X) 

a s r + a^ l + • • ■ + a r *• ' 

is a rational function in p and s with possible complex coefficients. 
For convenience, if we define 

6, = X*B t X 

w(p, s) = X*WX 

f(p, s) = b s r + 6,s' _1 + • • • + 6 r 

and as before 

g(p, s) = a s T + a^' 1 + • • • + a r 

equation (93) can be written as 

From the definition of a two variable reactance matrix, w(p, s) is a 
two variable positive function, and hence for any p with Re p > 0, 
w(p , s) is a positive function of s. 20 Consequently, for all s with Re s > 

Since equation (95) has to be satisfied for all s with Re s > and hence 
for arbitrarily small s, it can be seen from equation (93) that 

Re^^O 

for all p with Re p > 0. Hence, B r (p)/a r (p) is a positive real matrix 
and since W = — W f 

L«rJ LOrJ 

and thus B r /a r is a reactance matrix in p. 
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If instead of starting from the positive function f(p , s)/g(p , s), 
we start from d'f(p , s)/ds'/d v g(p , s)/ds", which is also a positive 
function 20 for all ^ v ^ r, the same arguments used in proving that 
B r /a t is a reactance matrix can be repeated to show that Bi/a t is a 
reactance matrix in p for all ^ i ^ r. 

Now to show that a,-+i/a, is a reactance function in p, we can use a 
similar proof based on the fact that d"~ l g{p , s)/ds'" 1 /d"g(po , s)/ds" 
is a positive function. 

Again, it can be seen from the fact that 

d'^gipo , s)/ds"- 1 /d'g(p , s)/ds' 
is a positive function. 20 that b 1 + 1 /6,- is a positive function satisfying 

M ■ -M' 

If X in equation (93) is chosen real 6, +1 /6< will be real for real p and 
hence XB i+1 X/XB { X for any real n X 1 matrix, X, is a reactance func- 
tion. 

To see that the zeros of a, are all simple and restricted to the imaginary 
axis: if any one of the a ; has a double zero on the jco axis or a zero off 
the ;'co axis, from the reactance nature of a i+1 /a { for all ^ i ^ r, 
all the a,- must have the same zero, and, consequently, W(p, s) will 
have an s-independent pole contradicting our original assumption that 
W has no such poles. 

We have thus proved all the claims of Theorem 3. 

appendix c 

Proof of the Miniviality of the Realization of W(p, s) 
in Both Variables 

In this appendix we show that the realization of W(p, s) that Sec- 
tion III gives is minimal in both the p and s variables. From the def- 
initions of S 8 [W(p, s)] and S p [W(p, s)], it can be shown that if 
W{p, s) is finite at p = <x> and s = oo, 

8.[W(p, s)] = 8.HP, s)] 

S p [W(p, s)] = 8MP, s)] 
where the two variable real polynomial 

n(p, s) = d (p)s k + d^s^ 1 + • • • + d k (p) (96) 
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is the least common denominator of all the minors of W(p, s). The 
form in which -q{p } s) is written in equation (96) immediately reveals 
that 



*.[W(p, •)] = k. 
And if Tj(p, s) is written as 

■nip, s) = c (p)p m + c,(p)p"- 1 + • • • + c m (p), (97) 

it can be seen that 

8,[W(p, «)] = m. 
c.i Minimum Elements 

Next we would like to find the minimum number of elements of each 
kind needed in the realization of W(jp, s). 

Lemma 1 states that k, the rank of T T -i{p), is equal to the s-degree 
of W(p, s), and the realization obtained there uses exactly k s-type 
elements. By equation (4) 

W(p, s) = Z tl (p) + 2 I2 (p)[2 22 (p) + *lJ- X Zu(p). 

Suppose that there exists a realization with k s-type elements, where 
k < k = rank r r _ x (p). Then, 

W(p, s) = z 11 (p) + 2 12 (p)[? 22 (p) + sl t .] _1 zj a (p) 

where the matrices z ia (p) and z 22 (p) are nXk and k X k , respectively. 
Then by equation (25), ZV^fo) = N(p)W(p) where 



#(p) = 



Sia(p) 
?ia(p)l22(p) 



is an nr X k matrix and hence, rank N(p) g k . Also, we have 

rank T r - X (p) ^ rank N(p) ^ ft < k = rank r r .,(p) 

which is a contradiction, and hence k = rank 7Y-i(p) = S s [W(p, s)] 
is the minimum number of s-type elements required in any realiza- 
tion. Now by repeating the same argument with a realization of 
W(p, s) where p-type elements are extracted instead of s-type ele- 



REACTANCE MATRICES 197 

ments, we can see that any realization must contain at least m p-type 
elements where m = 8 P [ (p, s) ] . 

c.2 Minimality o] the Realization in Section III 

We next discuss the minimality of the realization in both p-type 
and s-type elements. For the purpose of realization, the reactance 
matrix W{p, s) was decomposed as 

W(p, s) = z u (p) + z 12 (p)[z 22 (p) + shV'zlM (4) 

where 



Z(p) = 



(U) 



2n(p) Zi 2 (p) 

_-zM z 22 {p)_ 

can be realized as the impedance matrix of a lossless (n + k) port 
in the p-plane. W(p, s) is the impedance seen at the first n ports when 
the above (n + k) port network is terminated with unit s-plane in- 
ductors at its last fc-ports. Since k is the s-degree of W(p, s), the 
realization uses the minimum number of s-type elements. To show 
that the realization uses the minimum number of p-type elements, 
we have to show that 8[Z(p) ] = S p [W(p,s)]. For this we need a rela- 
tionship that exists between the least common denominator of the 
minors of W(p, s) and the determinant \z 22 (p) + sl fc |. 

Every minor of [z 22 {p) + si*] -1 can be expressed as (i(p, s)/ip{p, s) 
(See p. 21, of Ref. 12, Vol. 1) where /4p, s) and <p{p, s) are polynomials 
in s with coefficients from the field of rational functions in p. Further- 
more, 

<p{p, s) = | z 22 (p) + sl k | 

is a monic polynomial in s of degree k. 

Since W{p, s) has no p-independent or s-independent poles, every 
zero of tj(p, s) is a zero of ?(p, s) , and since k = S 8 [<p(p, s)] = 
8 8 [r)(p, s)], f(p, s) and rj(p, s)/d (p), which are monic polynomials in 
s with rational functions of p as coefficients, must be identical. Hence 

I «4> + .1* I - ^- W 

To show that 5[Z(p)] = 8 p [W(p, s)] (since we already know that 8[Z(p)] ^ 
8 p [W(p, s)]) it is sufficient to show that 8[Z(p)] ^ 8 v [W(p, «)]. To establish 
this inequality, consider the matrix S(p, s) defined by 

S(p, s) = [Zip) - sl n+k ][Z(p) + 8l n+k ]-\ (99) 
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When 8 = 1, S(p, s) is the scattering matrix of a lossless network, 
since Z(p) describes a lossless network and (see p. 184 of Ref. 11) 

8[Z(p)] = 8[S(p, 1)]. (100) 

Since S(p, 1) is para unitary (see p. 131 of Ref. 15) 

8[S(p, 1)] = 8[| S(p, 1) |]. (101) 

Equating the determinants of matrices on both sides of equation (99) 



l Qf m B \ I - I Z(p) - sl n+k | 



S(p, s) | = 



Using a formula from the theory of determinants (see p. 46 of Ref. 
12, Vol. I) 

(ffn sl n ) ~r Zi2\z 2 2 slk) Z12 | • | z 2 2 sl& | 

(Z n + Sin) + Z 12 (Z22 + Sl t ) _1 zl 2 | • | Z 22 + Sl k | ,■) 

= | TT(p, -s) - sl„ 1 | z 22 - sl k | 

I W( P ,S) + sl n | "|s 22 +Sl 4 f 

Now if \W(p, s) + sl„| is written as 

i "*■ s) + sl - i - rid < 103) 

where /i(p, s) is a real polynomial in p and s, since the left side of 
equation (102) is finite at p = oo 

9JLh(p, «)] ^ S P fo(p, s)] = 5 p [TF(p, ■)]. (104) 

Substituting equations (98) and (103) in equation (102), we have 

Hp, s) v(p, s) n(p, -s) d (p) 



yip, -s) h(p, s) do(p) v(p, s) 



(105) 



= h(p, -s) 
h(p,s) 

and by equations (100), (101), and (104) 

8[Z(p)] = 8[S(p, 1)] ^ 8 v [W(p, .)]. 

We have thus shown that 8[Z{p)] = 8 p [W(p, «)]. 

It should be noted that Z(p) is the impedance matrix of any loss- 
less coupling network in a realization of W(p, s), minimal in s, and 
hence we come to the important conclusion that if a realization of 
W(p, s) is minimal in one of the variables it is automatically minimal 
in the other variable. 
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